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The thermoelectric properties of cubic zincblend sihcon carbide nanowires (SiCNWs) with nitrogen impuri- 
ties and vacancies along [111] direction are theoretically studied by means of atomistic simulations. It is found 
that the thermoelectric figure of merit ZT of SiCNWs can be significantly enhanced by doping N impurities 
together with making Si vacancies. Aiming at obtaining a large ZT, we study possible energetically stable con- 
figurations, and disclose that, when N dopants locate at the center, a small number of Si vacancies at corners 
are most favored for n-type nanowires, while a large number of Si vacancies spreading into the flat edge sites 
are most favored for p-type nanowires. For the SiCNW with a diameter of 1.1 nm and a length of 4.6 nm, the 
ZT value for the n-type is shown capable of reaching 1.78 at 900K. The conditions to get higher ZT values for 
longer SiCNWs are also addressed. 
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I. INTRODUCTION 

Recent studies on the thermoelectric cooling and power 
generation by semiconductors as solid-state portable energy 
converters are prosperous, and there has been a great leap in 
both theoretical and experimental techniques^i^ Practical ap- 
plications require that the thermoelectric figure of merit ZT 
should be greater than 1.5 (e.g. Ref. 2). It is conceived 
that one of main approaches to promote ZT is to utilize low- 
dimensional semiconductor materialsi^ For one-dimensional 
(ID) quantum wires, as the scattering of confined electrons 
is dramatically avoided while the phonons are strongly scat- 
tered off the surface, the ZT could be significantly enhanced,'* 
which was supported experimentally in doped rough silicon 
nanowires (SiNWs),^ Bi2Te3 nanowires^ and Bi nanowires. ^•'^ 

On the other hand, silicon carbide (SiC) materials are at- 
tracting much attention for their thermoelectric properties. 
The nonmetallic semiconductor SiC shows excellent mechan- 
ical properties,^ chemical durability, and in particular, the high 
temperature stability. Owing to its wide bandgap and low in- 
trinsic carrier concentration, the semiconducting behavior of 
SiC can be kept at temperatures much higher than the case of 
Si, thus resulting in a higher operation temperature tolerance 
for SiC nanodevices.'*' Considering also the fact that an effec- 
tive recovery of waste heat from vehicle exhaust requires an 
operating temperature about 350 °C,^ and other applications 
such as high-temperature media-compatible flow sensors,— 
one can see that it is quite necessary to exploit the thermo- 
electric properties of low-dimensional SiC materials. 

As early as in 2003, Yoshida et al. have successfully doped 
N and B into 300 jim thick SiC film that contains Si and 
C vacancies, and observed that with N dopants, the power 
factor of the SiC film is raised nearly one order of magni- 
tude at 973 K.'2 For SiC nanowires (SiCNWs), the techniques 
of fabrication '° and measurements'-' have been remarkably 
improved recently, and theoretical analyses on structural 
properties,'"* mobility,'^ thermal conductivity,'^ etc., were 
also carried out for a few cases. These studies suggest that 
it is prospective to obtain nice thermoelectric performance on 
SiCNWs. To achieve this goal, the detailed simulations and 



studies on thermoelectric properties of SiCNWs are really es- 
sential, which is however sparse in literature. 

In the present work, we will focus on the 3C-SiCNWs 
doped with N impurities and vacancies along [111] direc- 
tion with an energetically favorable hexagonal cross section,— 
which are readily fabricated. '^''^^ We investigate the structures, 
electronic and phonon transmissions, and analyze the effect of 
defects on the thermoelectric performance of SiCNWs. It is 
found that the thermoelectric figure of merit ZT of SiCNWs 
can be significantly enhanced by doping N impurities together 
with making Si vacancies. An optimal doping strategy is sug- 
gested for both n-type and p-type SiCNWs. 

This paper is organized as follows. In Sec. II, the sim- 
ulation method and details will be described. In Sec. Ill, 
the structural stability of SiCNWs with defects will be dis- 
cussed. The electronic structures and transport properties of 
SiCNWs with various defects are shown in Sec. IV. In Sec. V, 
the phonon transport properties of SiCNWs will be analyzed. 
In Sec. VI, the effects of different defects on thermoelectric 
transport properties of SiCNWs are presented. Finally, a sum- 
mary will be given. 

II. CALCULATIONAL METHOD 

Structural optimizations are performed by the SIESTA 
code^i? which is based on the density functional theory 
(DFT)2° with norm-conserving pseudopotentials^' and linear 
combinations of atomic orbitals. The optimizations are spin- 
polarized, employing double-zeta polarized basis sets within 
generalized gradient approximation (GGA) expressed by PBE 
functional. The energy cutoff is 180 Ry, and the force tol- 
erance criterion is 0.04 tNIA for structural relaxation. The 
separation between neighboring nanowire surfaces is 15 ^. 

The electron transmission spectra are calculated by TRAN- 
SIESTA code,^^^ based on the Landauer-Biittiker and nonequi- 
librium Green's function (NEGF) formalism.ri2^ employing 
single-zeta basis sets with GGA and PBE. The energy cutoff 
is 100 Ry, and the convergence criterion of density matrix is 
0.005. The leads of the transport model are assumed to be 
pristine nanowires. To calculate the mean free path (MFP) 
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le and localization length ^ of single N dopant, the supercell 
between the leads is chosen to contain five unit cells with a 
total length of 38.5 A. The three middle units together are 
fully relaxed, where the single N dopant is in the central unit. 
Following the method introduced in Refs. |26 and T^, the 
average scattering resistance at different positions for single 
dopant is defined as 
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where Rs, Rc, Gi, Gq, pi and are scattering resistance, con- 
tact resistance, conductance with dopants, conductance with- 
out dopant, the weight of different doping positions, and the 
number of conducting channels, respectively. The mean re- 
sistance of wire and thus the MFP can be estimated linearly 
by 



Rc{E) 
(RsiE))' 



(3) 



where d is the average dopant-dopant separation estimated 
from a realistic doping density. Such a linear relation is valid 
when the wire length L is in the quasiballistic {L < 1^) and 
diffusive (l^ < L < ^) regimes. In the localization regime 
(L > ^), the resistance increases exponentially. ^''•^^ For 3C- 
SiC, the carrier concentration can be controlled in the range of 
10^^-10^^ cm"^ (Ref. |28]). In the present work we use d=40 
nm corresponding to a bulk doping density about 1.5x10 
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cm The localization length is calculated by^ 

m^l[NiE) + lUE). 



(4) 



To calculate the transport properties of the SiCNWs with 
different defect combinations of N dopants and vacancies, we 
suppose that the defects are distributed homogeneously in the 
units between the leads, where the length between the leads 
is taken as 3 units (about 2.3 nm) by default, except for the 
case of 6 units with particular specification. For phonon trans- 
mission, the force constant matrices are calculated by GULP 
code,^" which is based on the Tersoff model of empirical po- 
tential (TEP). To obtain dynamical matrices, the intra- 
atomic elements of the force constant matrices output from 
GULP are replaced by the data recalculated under the condi- 
tion of momentum conservation,-'"* where the atomic masses 
are then included. The precision criterion of the matrices is 
10^^. In terms of the dynamical matrices, the phonon trans- 
mission spectra are calculated within the Landauer-Biittiker 
and NEGF formalism similar to the case of electrons, follow- 
ing the method described in Refs. lf34l - [36ll . It is known that the 
Landauer-Biittiker formalism employs a finite central region 
for transport. In our calculations we use the nanowires with 
a fixed length to get the electronic and thermal conductance 
from the electron and phonon transmission spectra, where the 



length can be eliminated in calculations of the thermoelectric 
figure of merit ZT that can be obtained by ''''^'™ 
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where S is the Seebeck coefficient, T is the temperature, is 
the chemical potential, a^, and k^/j are the electronic con- 
ductance, thermal conductance of electrons and phonons, re- 
spectively, T(w) and T{E) are the phonon and electron trans- 
mission, respectively. S^cTe is the power factor Generally 
speaking, the closer of the chemical potential /i to the con- 
duction band minimum (CBM) or the valence band maximum 
(VBM) is, the higher carrier concentration of electrons or 
holes is.^^'^^ The nonlinear effects such as electron-electron, 
phonon-phonon and electron-phonon interactions are ignored 
for simplicity in the present worki^Sili 



III. STRUCTURAL STABILITY WITH DEFECTS 

In experimental and theoretical studies on SiCNWs, the 
most reported one is the cubic zincblend 3C(/3)-SiCNWs 
along [111] direction, because it can be deposited on 
Si I0ii449f4i, ^ typical model for simulation, we choose the 
ultrathin 3C-SiCNWs with a diameter of 1 . 1 nm and a hexag- 
onal cross section to investigate their structural and thermo- 
electric properties. As shown in Fig. 1(a), one unit cell of the 
pristine SiCNW contains 37 Si and 37 C atoms in pairs along 
the [111] direction, and 42 H atoms are covered on surface 
to saturate the dangling bonds of Si and C atoms. Follow- 
ing from the simulations on SiNWs^i^iiS it is necessary to 
include H atoms as surface passivation for electronic calcula- 
tions, while the H atoms could be omitted for phonon calcula- 
tions that produces a small deviation no larger than 3%i^^^ 

The properties of doping one N or B atom per unit cell into 
the 3C-SiCNW have been theoretically investigated in Ref. 
ll43ll . It reveals that it is more favorable to dope N atom sub- 
stituting C than substituting Si atom, the energetically most 
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FIG. 1. (Color online) The cross section of pure SiCNWH along the 
axial direction [ill], with the marks of different doping sites. The 
blue, dark red, orange and light green balls represent Si, C, H and N 
atoms, respectively. 

stable doping site of N is at the center of the SiCNW (site C), 
and the inner doping positions of B, C and D are energetically 
more favorable than the near-surface position E and surface 
positions A and F,"*^ which is in agreement with our results. 

We go on investigating the combined defects with N impu- 
rities and vacancies in the SiCNW. Such kinds of defects have 
been successfully implanted in SiC film,'^ but no details on 
the distribution of doping positions were given. Here, the sites 
of impurities in the SiCNW are marked in Fig. 1. For clarity, 
the case that dopes one N atom per unit cell substituting the 
C atom at the center (site C) is labeled as SiCNWNcH; the 
case that dopes two N atoms substituting C atoms at sites C 
and A is labeled as SiCNWNc^H; the case that dopes one N 
at site C and one Si vacancy at site A with saturated passiva- 
tion by H atoms is labeled as SiCNWNc V;|*H, where the case 
with unsaturated passivation leaving 3 dangling bonds per unit 
cell is labeled as SiCNWNcV;^*h(3), and leaving 1 danghng 
bond as SiCNWNcV;|^h(l); the case that dopes one N at C 
site and three Si vacancies at three corresponding A sites with 
saturated passivation is labeled as SiCNWNcVf^H; and so 
on. When vacancies are included, a supercell containing two 
units is employed where neighboring units have the interlaced 
vacancies at the corresponding symmetrical sites. Note that 
one dangling bond per unit cell in SiCNWNc^H is not left by 
the SiC, but by the N dopant at the surface site A. 

The formation energy E-^ averaged over one unit cell is cal- 
culated by^^ 

Ef{D) = E'°\D) - E'°\pure) - E.An.^i,, (11) 

where E*°*(D) and E*°*(pure) represent the total energy with 
and without the defects (D), respectively, An^ is the num- 
ber increment of atoms induced by the defects, /i^ is the 
chemical potential, and i = N,C, Si. At Si-rich, N-rich 
and H-rich limit, the chemical potentials are fisi = Ms""', 
Mc = MsiC ~ ^^Si, = and ^ih = mil'^^ re- 

spectively. The formation energies of several typical cases are 
shown in Table 1. 



TABLE I. The formation energies of SiCNWs with defects of N 
and vacancy under Si-, N- and H-rich conditions. 

NW with defect E-^ (eV) NW with defect E-^ (eV) 

SiCNWNcH -0.75 SiCNWNcVgH 1.86 

SICNWNcaH -1.38 SiCNWNcV|.*H -1.37 

SiCNWNcoH -1.36 SiCNWNcVf^H -2.51 

SiCNWNoVf h(l) -0.01 SiCNWNcVi'H -1.26 

SiCNWNoV5^h(3) 4.33 SiCNWNcVflH -1.45 

SiCNWNcVg.'h 3.89 SICNWHaVaH -0.86 



One may see that SiCNWNcV^*h(3), SiCNWNcVg.*h 
and SiCNWNcV^H have very large positive values of E^, 
implying that these structures with defects are unlikely to 
form because of too much formation energy needed4i^ 
SiCNWNcV;|*h(l) has an E-'^ very close to zero, showing 
that it is also unstable. All of the remaining cases have 
negative E^, which are energetically favorable and would 
be likely to be fabricated. The common ground of these 
cases is that all the Si and C atoms get saturated passiva- 
tion. It appears that the Si vacancies favor the surface sites 
where it is easy to passivate all the dangling bonds, thus en- 
abling those structures stable. Comparing SiCNWN^V^'H 
with SiCNWNcV^'H, in presence of the Si vacancy on the 
surface, the N dopant is obviously more favorable for the 
center site than the surface. Comparing SiCNWNc^H and 
SICNWNcdH, SiCNWNcVf'H and SiCNWNcV^^H, we 
note that in presence of one N dopant stable at the cen- 
ter, the E-^ differences between them are as small as 0.02 
eV and 0.11 eV, respectively, suggesting that their corre- 
sponding concentrations in the SiCNWs are close. How- 
ever, E-'^ of SiCNWNcVfpH is much lower than that of 
SiCNWNcVf^jH by 1.06 eV. It is clear that the location of 
Si vacancies prefers energetically the corner (site F) to the flat 
edge (site A). 

To estimate the MFP l^, and the localization length ^ for 
the case with single N dopant substituting one C atom, the 
electronic transmission spectra of the cases with one single N 
atom doped at different positions are calculated, as shown in 
Fig. 2(a). Except for the spectra very close to the CBM, the 
spectra associated with the three inner doping positions B, C 
and D are obviously higher than those of the surface and near- 
surface positions A, E and F. In Fig. 2(b), the MFPs for the 
structures with single N dopant of (i) the average of all six 
doping positions, (ii) the average of the positions B, C and D, 
and (iii) the position C are presented, where a dopant-dopant 
separation d is taken as 40 nm. The energetically stability of 
these three cases strengthens from (i) to (iii)'*-'. From Eqs. (3) 
and (4), we find that at E-EcBAf=0.12 with the conducting 
channels N=2, for case (i), (ii) and (iii), ^e=30.58, 64.78 and 
76.35 nm, ^=45.87, 97.17 and 114.53 nm, respectively; and 
at E-Eci3A/=0.27 with N=3, ie=19.84, 44.05 and 69.22 nm, 
C=39.69, 88.10 and 138.44 nm, respectively. It indicates that 
for single N dopant, the structures with energetically favor- 
able inner doping positions have the MFPs and localization 
lengths larger than those with the surface positions that are 
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FIG. 2. (Color online) (a) The electronic transmission spectra of 
SiCNWs with one single N dopant in the center of the supercell of 
5 units between the leads, doped at positions A, B, C, D, E and F, 
represented by dark yellow dash-dotted, green dotted, solid red, blue 
dashed, wine dash-dot-dotted and dark gray shot dashed curves, re- 
spectively. The black short dash-dotted curve represents the defect- 
free SiCNWH. (b) The MFP of single N dopant of three cases: (i) the 
average of all six doping positions (black short dash-dotted curve), 

(ii) the average of the positions B, C and D (blue dashed curve), and 

(iii) the position C (red solid curve), where d=40 nm. In both (a) and 
(b), the energy values of E=0.8 and 0.95 eV are marked as I and 2, 
respectively. 



FIG. 3. (Color online) The electronic band structure of (a) SiCNWH, 
(b) SiCNWNcH and (c) SiCNWNcVf'H (a supercell of 2 units). 
The zero point of energy is set to the VBM of SiCNWH. The corre- 
sponding Fermi levels are marked by dash-dotted lines. 
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IV. ELECTRONIC STRUCTURE AND TRANSPORT 

The electronic band structures of SiCNWH, SiCNWNcH 
and SiCNWNcVfH are shown in Fig. 3, where the zero 
point of energy is set to the VBM of the defect-free SiCNWH. 
From Fig. 3(a), one may see that SiCNWH has a direct gap, 
where the Fermi level locates in the band-gap. Fig. 3(b) man- 
ifests that the N dopant in SiCNW is an n-type defect, which 
significantly changes the bottom of the conduction band, and 
the Fermi level is cross this CBM4^ Therefore, the N impurity 
gives an essential contribution to the electronic transport. It is 
noted that the Fermi level of SiCNWNcH is 0.12 eV lower 
than the CRM of the defect-free SiCNWH, and the degener- 
acy of the bands of SiCNWNcH is also lower than that of SiC- 
NWH. Similar changes occur for the case of SiCNWNc Vf-'H 
shown in Fig. 3(c). The n-type defects introduce more elec- 
tron carriers and stronger scattering of electrons. 

Fig. 4(a) shows the density of states (DOS) of the defect- 
free SiCNWH, which is calculated with a spin-polarized code. 
It can be seen that the electronic band is spin unpolarized, and 



FIG. 4. (Color online) (a) The spin-polarized DOS of SiCNWH 
(black short dash-dotted curve) and the PDOS of Si (red solid curve) 
and C (blue dotted curve) atoms, (b) The spin-up DOS of SiCNWH, 
and the PDOS of one Si atom at positions A, B, C, D, E and F, re- 
spectively. The curve marks are in accordance with Fig. 2(a). The 
zero point of energy is set to the VBM of SiCNWH. 



the Si atoms contribute overwhelmingly to the DOS at the bot- 
tom of the conduction band. The peaks of the projected DOS 
(PDOS) of Si and C atoms totally overlap, and their relative 
magnitudes are inverse in the valence and conduction bands, 
indicating that the Si and C atoms bond with each other. The 
PDOS of one Si atom at different positions is presented in Fig. 
4(b). For the two peaks of DOS around 2.7 eV, the PDOS of 
the system with one Si atom at position B is almost coincident 
with that of Si at position C in the lower peak and with that 
of Si at position D in the higher peak. The PDOS of Si at po- 
sitions A and F in the higher peak are also nearly coincident. 
The PDOS of Si at inner positions B, C and D are significantly 
higher than those at outer positions E, F and A. 

It is interesting to note that for the peak of electronic trans- 
mission spectra from 0.8 to 1.2 eV in Fig. 2(a), the spectra of 
the system with N dopant substituting C atom at inner posi- 
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tions B, C and D are obviously higher than those at outer po- 
sitions E, F and A, which is quahtatively consistent with the 
PDOS of different Si atoms for the peak around 2.7 eV in Fig. 
4(b). When the N dopant substitutes one C atom in SiCNWs, 
the neighboring Si atoms connected by bonds will obtain ex- 
tra electrons from the dopants. The PDOS of the N dopant 
will spread to the valence and conduction bands similar to the 
PDOS of Si atoms. At the bottom of the conduction band, 
the Si atoms at inner positions contribute to the majority of the 
DOS, so the N defects at inner positions contribute electrons 
to the bottom of conduction bands more than those doped at 
outer positions. From these results one may deduce that for 
the peak of electronic transmission spectra at the bottom of 
conduction band, the spectra with N defect at inner positions 
could be higher than those at outer positions, which is just the 
results shown in Fig. 2(a). 

To probe the transport properties of SiCNWs with 
defects, besides the pure SiCNWH and SiCNWNcH 
that is energetically the most stable case for the Nc 
defect,43 SiCNWNcAH, SiCNWNcuH, SiCNWNcVf'H, 
SiCNWNcVf'H, SiCNWNcVf;^H and SiCNWNcVf],H 
are selected as relatively stable examples from the cases with 
different doping sites. 

The electronic transmissions of SiCNWs with various de- 
fects are given in Fig. 5. As shown in Fig. 5(a), for 
SiCNWNpH, at the bottom of conduction band the elec- 
tron transmission spectrum shrinks slightly from the spec- 
trum of SiCNWH, while at the top of valence band the hole 
transmission spectrum has a significant shrink from that of 
SiCNWH. For SICNWNcaH and SiCNWNcijH the trans- 
mission spectra go on shrinking from SiCNWNcH, and the 
shrink is extraordinarily evident for the hole transmission of 
SiCNWNcyiH that has one surface dangling bond left by the 
N dopant at site A. So the N impurities in the SiCNW induce 
stronger scattering of electrons, which overcomes the effect 
of the increase of electron carriers by the n-type doping, re- 
sulting in the decrease of electron/hole transmission. More 
dopants of N will lead to a greater decrease of transmission. 

The cases of SiCNWs with implantation of saturated Si va- 
cancies on surface in addition to the N dopant at center are 
shown in Fig. 5(b). The spectra of SiCNWNcVf^H and 
SiCNWNc Vgl^H with more Si vacancies are a bit lower than 
those of SiCNWNc Vf'H and SiCNWNcVf H owing to the 
stronger scattering by the vacancies. Comparing the spec- 
tra of SiCNWNcV;|*H, SiCNWNcVg^H, SiCNWNcVf'H 
and SiCNWNcVf),H in Fig. 5(b) with that of SiCNWNcH 
in Fig. 5(a), we find that in presence of the central N 
dopant, for the hole transmission with surface Si vacancies 
at either the corner (site F) or the flat edge (site A) site, 
and for the electronic transmission with surface Si vacan- 
cies at the corner F, the spectra change little from that of 
SiCNWNcH, revealing that the electron or hole scattering are 
not significantly strengthened by the vacancies. Consequently, 
the MFP and the localization length of SiCNWNcVf'H 
and SiCNWNcVfl^H would have values close to those of 
SiCNWNcH due to their similar transmission spectra. 

Fig. 5(c) compares the electronic transmission spectra of 
SiCNWNcVf'H with 3 and 6 units in length. The two 
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FIG. 5. (Color online) The electronic transmission spectra of SiC- 
NWs with various defects, where the black short dash-dotted curve 
always represent pure SiCNWH. (a) SiCNWNcH, SICNWNcaH 
and SiCNWNcuH are represented by red solid, green dot- 
ted and blue dashed curves, respectively; (b) SiCNWNcVi^H, 
SiCNWNcVfiH, SiCNWNcVf^'H and SiCNWNcVf^H are rep- 
resented by red solid, green dotted, blue dashed and dark yellow 
dash-dotted curves, respectively, (c) The electronic transmission 
spectra of SiCNWNc Vf.*H of 3 units (blue dashed curve) and 6 units 
(red dotted curve) for a comparison. 



spectra are quite close, and the spectrum of 6 units has a 
larger gradient at the edge of the bands. It reveals that for 
SiCNWNcV|.'H, the effects of increasing electron carriers 
and strengthening the electron scattering by the n-type defects 
get balanced so that the electronic transmission changes little, 
as the wire length increases. 

With the electronic transmission spectra, the electrical con- 
ductance (Te, the Seebeck coefficient S and the power factor 
S^(Te can be calculated by means of Eqs. (7) and (9). Gen- 
erally, the maximum of ZT is achieved when the chemical 
potential /i is around the CBM for electrons and the VBM for 
holes^^. The temperature dependences of ae and S are pre- 
sented in Figs. 6(a) and (b) for pure SiCNWH, SiCNWNcH, 
SiCNWNcV|*H and SiCNWNcVf|,H, respectively. These 
four cases are chosen owing to their relatively large transmis- 
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FIG. 6. (Color online) The temperature dependence of (a) the elec- 
trical conductance tJe and (b) the Seebeck coefficient S for pure 
SiCNWH (black short dash-dotted curve), SiCNWNcH (red solid 
curve), SiCNWNc-Vl'H (blue dashed curve) and SiCNWNcVf|,H 
(dark yellow dash-dotted curve), respectively. The chemical poten- 
tial fi is set to 0.675 eV, which is the bottom of the electron trans- 
mission spectrum of pure SiCNWH as indicated in Figs. 5(a) and 
(b). 



sion of electrons as seen from Figs. 5(a) and (b). Here fi is 
set to 0.675 eV, which is the bottom of the electron transmis- 
sion spectrum of pure SiCNWH. As shown in Figs. 6(a) and 
(b), the differences among SiCNWNcH, SiCNWNcVf^H 
and SiCNWNcVfl^H are small, unveiling that the effect of 
the increase of electron carriers and the enhancement of elec- 
tron scattering induced by more Si vacancies are nearly bal- 
anced. No matter the defects that include Si vacancies or not, 
CTe are significantly lower than those of pure SiCNWH due to 
the electron scattering by the N dopants. Moreover, it can be 
seen that for all four cases, as temperature increases from 200 
to 1000 K, the electrical conductance increases while the 
magnitude of negative Seebeck coefficient S decreases. When 
temperature is higher than 400 K, the changes of both and 
S become slow. 



V. PHONON TRANSPORT 

To investigate the properties of phonon transport. Fig. 7 
gives the phonon band structures of SiCNWH, SiCNWNcH 
and SiCNWNcVf-'H. It is uncovered that the degeneracy of 
phonon bands is lowered either by doping N at the center 
into SiCNWH, or by doping Si vacancies at the corner site 
F into SiCNWNcH, suggesting that the phonon scattering is 
strengthened. However, the number and the energy ranges 
of the phonon bands are nearly unchanged among these cases, 
which is in contrast to the defect bands in the gap of electronic 
band structures (Fig. 3). 

As shown in Fig. 8, the phonon transmission of pure SiC- 
NWH is equal to 4 at the low energy range close to zero, 
which corresponds to the four acoustic modes of phonons in 
nanowires. Compared with SiCNWH, the phonon transmis- 
sion spectrum of SiCNWNcVf.*H is obviously lower, and 
above 80 meV, the phonon transmission decreases near to 
zero. When the wire length increases from 3 to 6 units, the 
phonon transmission of SiCNWNcV|.*H is further lowered. 



(a) SiCNWH (b) SiCNWNcH (c) SiCNWN^VpH 

30 i 




FIG. 7. (Color online) The phonon band structure of (a) SiCNWH, 
(b) SiCNWNcH and (c) SiCNWNcVf'H, all of which are calcu- 
lated in a supercell of 2 units. 
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FIG. 8. (Color online) The phonon transmission spectra of pure SiC- 
NWH (black short dash-dotted curve), SiCNWNcVf'H with 3 units 
(blue dashed curve) and 6 units (red dotted curve), respectively. 
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FIG. 9. (Color online) The lattice thermal conductance Kph of SiC- 
NWs with various defects for temperature from to 1000 K. The 
symbols of different curves are totally the same as Figs. 5(a) and (b). 



In contrast to the results in Fig. 5(c) that the electronic trans- 
mission of SiCNWNcV|.*H change little from 3 to 6 units 
in length, there are no factors to balance the strengthening of 
phonon scattering by the defects, and the phonon transmission 
of SiCNWNcV|,*H decreases significantly as the wire length 
increases. 

Fig. 9 presents the lattice thermal conductance Kph for 
various cases in the temperature range of 0-1000 K. Kph 
first goes up rapidly and then increases slowly when tem- 
perature is higher than 300 K. It is seen that the Kph de- 
creases as the defects of N or Si vacancy increase. The dif- 
ferences of K^h between SICNWNcaH and SiCNWNcoH, 
SiCNWNcVf H and SiCNWNcVf'H, SiCNWNcVf^H 
and SiCNWNcVf|,H are small, which are 0.01 nW/K, 0.03 
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nW/K and 0.01 nW/K at 1000 K, respectively. It indicates that 
the same doping concentration at different doping sites leads 
to very similar thermal conductance. Comparing Fig. 6 with 
the corresponding cases in Fig. 9, one may see that in the 
presence of N dopants or Si vacancies, both de and Kph will 
decrease from the values of pure SiCNWH, but the amplitude 
of Kpfi is decreased more remarkably than that of a^. 



VI. THERMOELECTRIC TRANSPORT 

Combining the results of electron and phonon transport, the 
thermoelectric figure of merit ZT is calculated as a function of 
the chemical potential by Eqs. (5)-(10) at 300 K and 900 
K, respectively, as shown in Fig. 10. Comparing Figs. 10 and 
5, the maximum of ZT occurs when the chemical potential 
reaches the VBM for p-type SiCNWs or the CBM for n-type 
SiCNWs, where the electronic transmission changes dramati- 
cally. 

For n-type SiCNWs, at 300K or 900K, it is clear that in 
Figs. 10(a), 10(b), 10(d) and 10(e), the cases of SiCNWH^H, 



SiCNWNcV|,'H and SiCNWNcVf>H can induce obvi 



pn and SiCNWNcVf>: 
ous enhancements of the ZT. For other cases the ZT is 
even lower than that of SiCNWH. Among SiCNWNcH, 
SiCNWNcVf'H and SiCNWNcVf|,H, their relative magni- 
tudes of ZT change as the length increases from 3 to 6 units, 
as shown in Figs. 10(c) and 10(f). For a length of 6 units, the 
n-type ZT maximum of SiCNWNcH, SiCNWNcVf^'H and 
SiCNWNcVf^H are 0.45, 0.55 and 0.43 at 300 K, 1.37, 1.78 
and 1.67 at 900 K, respectively. Thus, the ZT maximum for 
n-type is given by SiCNWNc Vf.*H, containing the defects of 
one N dopant at the center and one Si vacancy at the corner 
per unit cell. Neither more N dopants such as SiCNWNp dH 
nor more Si vacancies such as SiCNWNcVf|,H are benefi- 
cial. In general, to avoid the ZT dropping from the maximum 
of SiCNWNcV|^*H to the value of SiCNWNcVf|,H or oth- 
ers, it is preferable to dope N impurities limited at the center 
and a small quantity of Si vacancies limited at the corners (F 
sites), which is energetically favorable. Such a distribution of 
Si vacancies at the corners strengthens the phonon scattering 
dramatically while it influences little on the electron scatter- 
ing. 

For SiCNWNcH and SiCNWNcVf^'H, it was shown in 
Fig. 3 that their Fermi levels are about 0. 1 eV lower than the 
CBM of the defect-free SiCNWH. The Fermi level equals to 
the chemical potential at zero temperature. As manifested in 
Fig. 10, because the highest peaks of ZT for n-type accumu- 
late around the CBM of SiCNWH, the corresponding chemi- 
cal potential at finite temperature will lead to a ZT value very 
close to the maximum of ZT for n-type. 

For p-type SiCNWs, among the cases in Figs. 10(a), 10(b), 
10(d) and 10(e), the maximum of ZT of SiCNWNcpH is 
larger than that of SiCNWNc^H. Considering the tiny dif- 
ference of formation energy E-' between SiCNWNcyiH and 
SICNWNcdH, the larger ZT of SiCNWNccH is hard to 
maintain. On the other hand, the differences of ZT among 
SiCNWNcV^'H, SiCNWiVcVfjiH, SiCNWNcVf*H and 
SiCNWNc Vfl^H are quite small, while the energetically most 
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FIG. 10. (Color online) The thermoelectric figure of merit ZT of 
SiCNWs with various defects as a function of chemical potential fi 
at (a)-(c) 300 K and (d)-(f) 900 K. The length of SiCNWs in (c) and 
(f) is 4.6 nm of 6 unit cells, and the length of other cases is 2.3 nm 
of 3 unit cells. The symbols of different curves are totally the same 
as in Figs. 5(a) and 5(b). On both sides of the zero point of chemical 
potential, the ZT value of the first peak of each curve is listed in the 
bracket after the figure legend of the corresponding curve. The left 
(right) value in the bracket corresponds to the first peak of p-type 
(n-type) on the left (right) side. 
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favorable one is SiCNWNcVf|,H. All these four cases have 
significant ZT enhancements in Figs. 10(c) and 10(f), and rel- 
atively, the cases with more Si vacancies have larger ZT. For 
the best two cases, SiCNWNcVfj^H and SiCNWNcVf|,H 
with a length of 6 units, the p-type ZT maximum is 0.47 and 
0.46 at 300 K, 1.66 and 1.49 at 900 K, respectively, which 
are obviously larger than the values of SiCNWH, as shown in 
Figs. 10(c) and 10(f). Therefore, it is preferable to dope N im- 
purities stable at the center and a large quantity of Si vacancies 
not only at the corners (F sites) but also spreading to the flat 
edges (A sites). However, the p-type performance depends on 
the ingredient of the p-type dopants which is not addressed in 
the present work. 

It should be noted that the selected cases cannot cover all 
positions and concentrations of the defects. However, the 
performances affected by various defects are reflected in the 
ZT curves of these selected interesting cases as functions of 
/i. For instance, as displayed in Figs. 10(a) and 10(d), for 
the n-type the peaks of ZT for SiCNWNcH, SiCNWNc^H 
and SICNWNcdH appear at ^ = 0.675, 0.695 and 0.725 
eV at 300 K, 0.595, 0.635 and 0.665 eV at 900 K, re- 
spectively. Compared with SiCNWNcH, SiCNWNc^H and 
SiCNWNc£)H have more N dopants leading to a higher car- 
rier concentration of electrons and a larger chemical poten- 
tial. The ZT maximums are listed in the brackets in the fig- 
ure legends. It can be seen that for the n-type SiCNWNc-H, 
as the chemical potential /i increases, the ZT values on the 
right side of the highest peak are close to the ZT maximums 
of SiCNWNc^H and SiCNWNcnH. Thus, as seen in Figs. 
10(c) and 10(f), even if the ZT curves of SiCNWNc^H and 
SiCNWNc £)H are not shown, for the cases with different dop- 
ing concentrations of N, the corresponding ZT could be ap- 
proximately estimated from the ZT curve of SiCNWNc-H. 

From Figs. 10(c) and 10(f), it is clear that the ZT of 6 units 
is larger than the ZT of 3 units for all the cases under inter- 
est. As an example, the electronic and phonon transmission 
spectra of SiCNWNcV|.*H with different lengths have been 
illustrated in Figs. 5(c) and 8. From the discussions in Figs. 3, 
5(c), 6, 8 and 9, we know that for the cases of n-type doping, 
the effect of strengthening the electron scattering is balanced 
by the increase of the electron carriers for electronic trans- 
mission, while the strengthening of the phonon scattering is 
unbalanced. Thus, as the wire length increases from 2.3 to 4.6 
nm, the electronic transmission change little while the phonon 
transmission decreases significantly, leading to the increase of 
ZT value. 

In Sec. Ill, the localization length of the case with N dopant 
was estimated to be greater than 100 nm at the bottom of con- 
duction band. Doping Si vacancies at the corners into the 
N doped SiCNW will not change the localization length ob- 
viously, for SiCNWNcV|.*H has an electronic transmission 
spectrum close to SiCNWNcH. With a wire length smaller 
than the localization length, the electronic resistance will in- 



crease along with the wire length linearly. • Moreover, it 
has been reported in the case of SiNWs with a diameter of 1 .2 
nm that the thermal conductance keeps decreasing dramati- 
cally before the length increases over 50 nm, and after this 
length the thermal conductance decreases almost linearly.— 
In this sense, for the nanowires such as SiCNWNcV|.*H, the 
tendency of ZT enhancement is probable to be kept not only 
for the wire length from 2.3 nm to 4.6 nm, but also for a wire 
length of about 50 nm. As a result, there is still a large room 
to expect a further enhancement of ZT when the length in- 
creases from 4.6 nm, and a large ZT can be maintained within 
the localization length. 

VII. SUMMARY 

In summary, the structural stability and the thermoelectric 
properties of SiCNWs with N dopants and vacancies are in- 
vestigated by means of the density functional calculations. 
Our studies show that the central (C) site for N dopants is 
energetically the most stable. When the defects contain both 
N impurities and vacancies, the most favorable configuration 
is that the N impurities are at the center while the Si vacancies 
locate at the corners (F) with all the dangling bonds passi- 
vated. Aiming at obtaining a large thermoelectric figure of 
merit ZT, with N dopants at the center, a small quantity of 
Si vacancies limited at the corners is most favored for n-type, 
while a large quantity of Si vacancies spreading to the flat 
edge (A) sites is most favored for p-type wires. 

For the SiCNW as quantum wires along [111] with a diam- 
eter of 1 . 1 nm and a length of 4.6 nm, the ZT maxima at 900 K 
are found to be about 1.78 for n-type and 1.66 for p-type, re- 
spectively. The ZT maximum of n-type can be reached by the 
present n-type doping of N impurities and Si vacancies, but 
for p-type it depends on whether or not there are proper p-type 
dopants. As long as the decrease of phonon transmission with 
the increase of length is far more rapid than the decrease of 
electronic transmission, higher ZT values may be expected for 
SiCNWs longer than 4.6 nm. Based on the reduced model of 
quantum wires for simulation, the present findings may shed 
light on the enhancement of the ZT values by taking advan- 
tage of N dopants and surface Si vacancies in SiCNWs, that 
is especially favorable for the thermoelectric performance at 
high temperature in applications. 
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